Phenomenology of Conduction in Incoherent Layered Crystals 
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A novel phenomenological approach to the analysis of the conductivities of incoherent layered 
crystals is presented. It is based on the fundamental relationship between the resistive anisotropy 
o 'abl^c and the ratio of the phase coherence lengths in the respective directions. We explore the 
model-independent consequences of a general assumption that the out-of-plane phase coherence 
length of single electrons is a short fixed distance of the order of interlayer spacing. Several topics 
are discussed: application of the scaling theory, magnetoresistivity, the effects of substitutions and 
the intermediate regime of conduction when both coherence lengths change with temperature, but 
at different rate. 
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I. INTRODUCTION 

The normal state conductivity of highly anisotropic 
layered crystals such as underdoped high-T c cuprates ex- 
hibits a number of unusual features widely discussed in 
literature. Metallic in-plane (a a ) and non-metallic out- 
of-plane (<t c ) DC conductivities coexist in a broad ranges 
of temperature and doping 1]. A prominent feature of 
underdoped layered cuprates is strongly temperature de- 
pendent anisotropy a a /a c , which in many cases shows no 
tendency to saturation even at low temperatures. Op- 
tical conductivity a c (u) is approximately frequency in- 
dependent over a wide range of frequencies. These and 
other findings have led to conclusion that the out-of-plane 
transport in these crystals is incoherent 0, IE 0, IE IE ■ 

The microscopic models based on the assumption of in- 
coherence of the c-axis transport such as that in Rcf. 8] 
and othersplQ corroborate this conjecture, because the 
properties of a c obtained within these models reproduce 
some of the features of the real systems. Presently, how- 
ever, the microscopic models are still too restrictive and 
idealized to provide a comprehensive description of the 
phenomenon of incoherent conduction and to be a versa- 
tile framework for the analysis of the experimental data. 

It would be highly useful to have a phenomenology of 
incoherent transport that would rely minimally on the 
microscopic models and was based mainly on the ex- 
ploration of the consequences of a general assumption 
that the out-of-plane phase coherence length of the charge 
carriers is a short fixed distance of the order of inter- 
atomic distances. A novel approach presented here is 
based in large part on the relationship between resis- 
tive anisotropy a a /a c and phase coherence lengths of the 
charge carriers. 

The scope of the paper can be summarized as follows: 
In Sees. II and III we have derived the relationship be- 
tween the anisotropy and phase coherence lengths and 
compare the results with the solvable microscopic mod- 
els. In Sec. IV the scaling theory is applied in order to 
develop a more comprehensive description of conductiv- 
ities than that offered by existing microscopic models. 



When the out-of-plane phase coherence length £ VyC is 
fixed, the resistive anisotropy is T-dependent and reflects 
the T-dependence of the in- plane coherence length l Vta , 
namely r\ = (aa/ac) 1 / 2 cx £ Vta . Thus, the analysis of the 
dependence of conductivity on the anisotropy, a vs 77, 
as opposed to conventional — conductivity vs. tempera- 
ture — approach, gives an opportunity to gain an insight 
into the dependence of the conductivities on the coher- 
ence length. An extention of the one-parameter scal- 
ing hypothesis |E E3 on the conductance of the phase- 
coherent volume leads to an experimentally verifiable 
prediction that there may exist a unifying description 
of conductivity of incoherent crystals at different doping 
levels of the form a /a — f(T]/fj), where f(y) is a uni- 
versal function for a given class of crystals, and a and f\ 
are doping-dependent normalization constants. The scal- 
ing theory allows to obtain the functional form of this 
dependence for insulating and metallic branches of the 
trajectories. For the metallic branch a a b ~ ln^ and for 
the insulating branch p ao ~ blip. At low temperatures 
this translates into logarithmic temperature dependence 
of resistivity observed earlier The logarithmic 

dependence of conductivity is also shown to be present in 
Y\^ x Pr x Ba2Cu'iO x crystals. 

In Sec. V we discuss magnetoresistivity. An impor- 
tant result is that in the temperature range where one 
can observe the magnetoresistance caused by quantum 
interference 13], the relationship between anisotropy of 
incoherent crystals and the in-plane coherence length: 
a a b/a c oc I ah , can be verified experimentally. Since the 
crossover field of magnetoresistance B v also depends on 
the phase coherence length, B v oc ab , a combination 
of seemingly unrelated quantities - B v a ao /a c - should re- 
main constant, even though each of the factors, B v and 
anisotropy, strongly change with temperature. 

Section VI is devoted to the effects of elemental sub- 
stitutions and disorder in incoherent crystals. The treat- 
ment addresses sometimes puzzling effect of substitu- 
tions for Cu and radiation damage on resistivities. In 
Sec. VII we discuss the " semicoherent crystals" in which 
both coherence lengths are changing with temperature, 
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but at different rates, so that the anisotropy is still 
T— dependent, but differently than in fully incoherent 
crystals where l V:C = const. Section VIII offers some sug- 
gestions on the type of future experiments that can shed 
light on the nature of confinement in cuprates. 



II. THE RELATIONSHIP BETWEEN 
ANISOTROPY AND PHASE COHERENCE 
LENGTH 

Kubo-Greenwood formula relates the diagonal compo- 
nents of the conductivity tensor to the matrix elements of 
the respective components of the position operator |l4j: 



if £2 i^PCfo 



f a )u>6(e a - 6/3 - hw). 

Here {x, y, z} are the principal axes of the crystal, O is 
the volume of the system, oj is the frequency, and the 
static conductivity a xx = lim w _»o <j xx {oj). x a p are the 
matrix elements of the x— component of the position op- 
erator between states of energy e a ,ep, whose probabilities 
of occupation in thermal equilibrium are given by f a , fp- 
Two other diagonal components of the conductivity ten- 
sor, o yy and a zz , are determined, respectively, by the 
matrix elements y a f3 and z a p. 

The directional dependence of the conductivities is de- 
termined entirely by the matrix elements of the position 
operator. Therefore, it is obvious that one can always 
find the length scales {A^, X y ,X z } such that 

Mx/\ x \(3)\ 2 = E iHy/Wi 2 = £iwa*i/3>i 2 . 



OC0 



a/3 



a/3 



The summation here is over the states relevant to con- 
ductivity (1). Thus, the anisotropy is determined by the 
ratio of these length scales: 



O ~xx 
a yy 



X y 



®~xx 
0~zz 



X x 



(2) 



What is the physical nature of these characteristic 
lengths? If we look at the ways the Eq. (1) can be re- 
duced to the standard result of kinetic theory [lHIi5| . the 
matrix elements in (1) are obtained by integrating over 
the phase-coherent volume ( which is, roughly speaking, 
a rectangular block with the sides equal to the phase co- 
herence lengths in the respective directions). Therefore, 
at least in the case of Fermi liquid type system, the length 
scales determining anisotropy in Eq. (2) are the phase 
coherence lengths. 

A deeper insight into the meaning of the length scales 
in Eq.(2) can be gained from Thouless' concept of con- 
ductance of a microscopic block. The conductance of a 
block whose edges are along the pri ncipal axes of the 
conductivity tensor is given by 

e 2 dN 



ne principal s 



where i = {x, y, z}, dN/dE is the total number of states 
inside the block per unit energy and (AE)i is the mean 
fluctuation in energy levels caused by replacing periodic 
by antiperiodic boundary conditions normal to the direc- 
tion of the current. 

The sensitivity of the energy spectrum to the bound- 
ary conditions depends on how well the wave function 
retains its phase coherence along the path between the 
boundaries. We can choose the sides of the block such 
that, on average, the random phase acquired due to an 
inelastic interaction along the way between the two op- 
posite boundaries is the same, of the order of tt, for all 
three pairs of the block boundaries. Then, the effect of 
the imposed phase difference between the boundaries on 
the energy spectrum of such a block is isotropic; namely, 
(AE) X = (AE) y = {AE) Z . According to Eq. (3) the 
conductance of such a block is isotropic. This choice of 
the sides of the block corresponds to the definition of 
the phase coherence length £ v i : the distance over which 
electrons lose phase coherence [l^. Thus, the conduc- 
tance g Vt i of the phase-coherent volume of an anisotropic 
medium is isotropic: 



(4) 



h dE 



(AE) % 



The conductances of neighboring phase-coherent vol- 
umes add up according to Ohm's law ^3|. Therefore, 
a macroscopic block {L x , L y , L z } obtained by fitting to- 
gether N 3 phase-coherent volumes (L x /£ v . x = L y /i v ^ y — 
Lz/i^.z = N 3> 1) also has an isotropic conductance 
G « Ng v , which can be expressed in terms of the com- 
ponents of the conductivity tensor: 

Go~ xx L y L z o~yyL x L z o~ zz L y L x . . 

= — r 1 — = —j — = — j*—. (5) 

Lix i-"y i~i z 

This leads to the following relationship between conduc- 
tivities: 

0~xx <P,z . 0~xx <p,x 

a yy ^ip,y a zz l ip,z 

The Thouless' definition of conductance [Eq. (3)] is a 
very general approach to dissipation, because it relates 
conductance only to the statistical properties of the en- 
ergy spectrum. Therefore, the conjecture expressed by 
Eqs. (4) and (6) may have a general applicability, re- 
gardless of the type of conduction. 



A. Comparison with Fermi liquid 

A characteristic property of conventional Fermi liquids 
is that the phase coherence of the wave function extends 
over many unit cells in all directions. The coherence 
lengths in different directions change with temperature 
at the same rate, determined by the inelastic relaxation 
time which is equal, or scales with, the decoherence time 
t„. As the result, the anisotropy of Fermi liquids is tern- 



(3) perature independent 
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Specificaly, in a Fermi liquid with ballistic transport 
, t i oc VpjTu,, so that Eq. (6) reduces to: 



(7) 



Here vp,i is an average component of the Fermi veloc- 
ity and Tip is assumed to be independent of the veloc- 
ity. Equation (7) coincides with the well known result 
for Fermi liquids which follows from the solution of the 
kinetic equation under the same assumption, that the de- 
coherence time is independent of the velocity. In the case 
of the parabolic dependence of the energy upon momen- 
tum, Eq.(7) reduces to the ratio of the effective masses: 



(8) 



Thus, Eq.(6) is consistent with the Fermi liquid result 
for anisotropy. At first glance this statement appears 
counter-intuitive since the conductivity of Fermi liquid is 
proportional to the first power of inelastic relaxation time 
and, consequently, the first power of £ v . But a closer look 
at the expression for anisotropic conductivity ^jj imme- 
diately clarifies this misapprehension: an oc e (vp^T^ 
can be written as an oc e t%,i/T v . Now au directional 
dependence of the conductivity is incorporated into lr f . 
The relaxation time in the denominator cancels out in 
the anisotropy, which leads to Eq. (6). 

The comparison with Fermi liquids also allows us to 
clarify the definition of coherence lengths in Eq. (6). It 
should be understood that 0? { in Eq. (6) is the mean 
square average of the respective coherence length. 

In case of diffusive motion oc DiT v ) Eq. (6) re- 
duces to the Einstein relation <Ji/o~j — Di/Dj. The re- 
verse procedure of obtaining Eq. (6) from the Einstein 
relation is not as straightforward and requires a micro- 
scopic model of conduction. See an example in Ref . [T3| . 



III. 



INCOHERENT CRYSTALS 



Relationship (6) provides an effective tool for develop- 
ing a phenomenological description of electric transport 
in incoherent layered crystals. The nature of interlayer 
decoherence does not concern us here. It may be the 
result of the intralayer processes [when the rate of tran- 
sitions between layers u± is small in comparison with 
the decoherence time (u>j_t v < 1)] or interlayer transi- 
tion itself (if a random phase of the order of it is acquired 
during the transition from one plane to another). In ei- 
ther case, the essential point is that the phase-coherent 
volume of such a crystal contains only one layer (bilayer 
in the case of cuprates like YBa-iCu^Ox). Correspond- 
ingly, the out-of-plane phase coherence length, hereafter 
denoted as £o, is the T— independent distance of the or- 
der of interatomic (interlayer) spacing. Obviously, the 
coherence length cannot be arbitrarily small: it cannot 



be smaller than the size of the atomic orbitals. The im- 
mediate consequence of Eq. (6) is that in incoherent 
crystals 



Cab 

<Tc. 



f 2 
f 2 ' 



(9) 



Here, for simplicity, we consider isotropic planes a x 



'vv 



= a ab , and i, 



<-v,y 



Since £ v monotonically increases with decreasing tem- 
perature, so does the anisotropy. Thus, a fundamental 
feature of incoherent crystals is the anisotropy which re- 
flects the T-dependence of the in-plane phase coherence 
length. A consequence of Eq. (9) is that the out-of- 
plane normal state conductivity of incoherent crystals is 
completely determined by the in-plane conductivity and 
phase coherence length: 



I 2 ' 



(10) 



An important theoretical result that confirms Eqs. (9) 
and (10) is the work of Graf, Rainer and Sauls They 
modeled a layered metal as a stack of 2D conducting 
planes coupled via incoherent interplane scattering of 
charge carriers. The planes were treated as 2D Fermi 
liquid. In the limit when the decoherence time t v is 
the same for in-plane and out-of-plane conductivities, the 
anisotropy obtained inQ is given by: 



Q~ab 



Ad 2 • 



(11) 



where d is the interplane distance and vp is the Fermi 
velocity of the circular Fermi surface. If we define tz, in 
Eq. (9) as the mean square average over the Fermi circle: 



c = v 



F'cp 



tt/2 



tt/2 



2tt 



■ cos 



Eq. (11) reduces to Eq. (9) with £q = d. If we take the 
planes to be conventional two dimensional Fermi liquids 
with a ab = qv F T m ~ qvpT^ = q'£ 9 , Eqs. (10) and (11) 
give a c = q'£l/£ Vl or a c a a b = const = q a l\. 

The comparison of Eq. (9) with Eq. (11) seems to in- 
dicate that when we apply Eq. (9) to the analysis of the 
experimental data in Y BaiCu?,O x , for example, we have 
to take £o equal to the total distance between the bilay- 
ers {12 A). However, the real crystals are not necessarily 
can be adequately described by a model of stacked metal 
planes. The analysis of magnetoresistivity given in Sec- 
tion VI indicates that £q in YBaiCuj,O x may be equal 
to the half of the distance between neighboring CuO^ 
bilayers (6 A). The decoherence, apparently, is taking 
place over the distance between a CuOi bilayer and the 
neighboring CuO chain layer. The exact value of £q will 
become clear only when we understand the mechanism 
of decoherence in cuprates. More details on this matter 
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are given in Section VII. Hereafter I will take £q equal 
to the half distance between the neighboring CuOi bi- 
layers. However, most of the results presented below do 
not depend on the exact value of £q. 

Below we consider experimental data that corroborate 
the validity of Eqs. (9) and (10). It concerns two types 
of incoherent crystals: insulating and optimally doped 
cuprates. 



A. Hopping conduction in insulating layered 
crystals 

The first example which allows to test the validity of 
the conjecture (6) and its consequences (9) and (10) is the 
anisotropic hopping conduction. The coherence length of 
localized carriers is equal to the hopping distance. If the 
hopping distances in different directions are temperature 
independent, as is the case in the nearest neighbor hop- 
ping regime, or if they change with temperature at the 
same rate, which corresponds to anisotropic 3D variable 
range hopping, the anisotropy remains T— independent 
even though it may be very large. Exponentially strong 
T— dependence of the conductivities cancels out. It 
should be noted that the result equivalent to Eq. (6) was 
obtained earlier for the critical network model of the hop- 
ping transport [l^. It was shown that the anisotropy is 
given by the square of the ratio of the correlation lengths 
of the critical network. 

In incoherent layered insulating crystals like 
PrBa2Cu^O'j-s the localized states are two-dimensional, 
i.e. comprised of the orbitals that all belong to one 
bilayer 18J. In the variable range hopping (VRH) regime, 
the in-plane hopping distance increases with decreasing 
temperature, while in the out-of-plane direction the 
carriers advance in fixed steps equal to the distance 
between the neighboring bilayers. Thus, the resistive 
anisotropy of such crystals increases at low temperatures 
when VRH sets in. To determine the T— dependence of 
the anisotropy, let us consider the conventional picture 
of VRH where the average in-plane hopping distance 
can be found by maximizing the hopping probability 



P(R) cx exp <{ -2- - 



A 



X Af(R 2 - B%)T 



(12) 



Here A is the localization length, JV = const is the 2D 
density of states, and A a numerical coefficient. The only 
modification of the traditional treatment of VRH is the 
denominator R 2 — i?, 2 , instead of R 2 . This takes into ac- 
count that two localized states with close energies cannot 
be found closer to each other than a certain distance Rq of 
the order of the localization length A. Indeed, if the two 
states with close energies strongly overlap, the phonon 
interaction that causes hopping will also hybridize them 
and push apart the energies of the new states. Equation 
(12) assures that in the limit of high temperatures the 
hopping distance does not decrease below the value R 
which corresponds to the nearest-neighbor hopping. In 



the limit T — > 0, Eq. (12) crosses over into conventional 
Mott form. Strictly speaking, Eq. (12) is valid only for 
R 3> Rq ~ 2A which corresponds to VRH regime. In 
this limit the average hopping distance £h, which in the 
hopping regime is equal to the in-plane coherence length, 
is determined by the maximum of hopping probability 
(12): 



th~R 



2< 
3i? ' 



R 



(—) 

\J\fT J 



1/3 



Then, according to Eq. (9), 



a ab £\ ^ m 2 ( XA \ 



2/3 



( 2 



3£ 2 



(13) 



This T-dependence of the anisotropy (a + 6T~ 2 / 3 ) was 
observed in insulating PrBa2Cu^Oj-s and strongly un- 
derdoped superconducting Yj -^Pr^Bar>Cu^Q7-x\l8^. In 
PrBa 2 Cu 3 7 ~5, a w 123 and b « 171 K 2 I Z . It is in- 
structive to show that Eq. (13) not only determines cor- 
rectly the temperature dependence of the anisotropy, but 
is quantitatively correct as well. The exponential temper- 
ature dependence of conductivity is determined by the 
maximum of the probability (12) (a ~ P m ax) at R w R: 



3R 

a ~ exp \ — — f = exp 



Tm 
T 



1/3- 



;T M = 



27 A 



AW 
(14) 

where Tm is the Mott parameter. The anisotropy (13) 
can be presented as 



Q~ab 



Thus, 



T, 



2/3 



27£| 
A 3 



T a = 



XA 

w 



T„ 



(15) 



(16) 



In Ref.0 we reported that T M « 300 K and T a = 
fe 3 / 2 w 2236 K. This gives A = 3£ (T Q /T M ) 1/3 « 6£ - 
The constant term a determines the cutoff length Rq — 
4(3a/4) 1 / 2 w 9.64 = 1-6A. This value of the cutoff R 
is within the expected range A < Rq < 2A. For £ w 6A, 
we get A « 36A 



B. Resistivity of optimally doped high-T c cuprates 

As I already mentioned above, the assumptions that 
the in-plane conductivity of incoherent layered crystal is 
the Fermi liquid type: a a b ~ e 2 nTi n /m and the coherence 
length is proportional to the inelastic relaxation time Tj n 
lead to conclusion that the product <J c (J a b must be tem- 
perature independent. Even though there is a great deal 
of evidence that optimally doped cuprates are incoherent, 
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the temperature dependence of in- and out-of plane con- 
ductivities do not follow this rule. Since the anisotropy 
of such crystals directly determines the in-plane phase 
coherence length, Eq. (9), one can obtain the depen- 
dence conductivity vs. coherence length, a a b(£ v ), sim- 
ply by plotting the conductivity data against anisotropy. 
In Ref. 20] the in-plane conductivity of nearly optimally 
doped Y Ba2Cu30e,g3 crystal was analyzed in this way 
over the temperature range 90 < T < 300 K. The critical 
temperature of that sample was slightly below 90 K . The 
in-plane conductivity can be very accurately described 
as: 

ffai = g(^-0; ( 17 ) 

Here £ v = £ (p c / Pab) 1/2 and £ « 2lA. This linear de- 
pendence da a b/d£ ip = const indicates ballistic motion of 
the quasiparticles. Indeed, the quasiclassical conductiv- 
ity (Jab oc t (r is the relaxation time of the distribu- 
tion function). For ballistic motion d£ v /dT ip = vp (t v is 
the decoherence time) and if r oc t v , then dcr a b/d£ v oc 
da a b/dr v — const. According to Eq. (10), the correspond- 
ing a c is then given by 




Note that da c /d£ v > for £ lf> < 2£, and da c ld£ v < 
for l v > 2£. Since the phase coherence length mono- 
tonically increases with decreasing temperature, a c is 
metallic at high temperature (£ v < 2£), reaches a max- 
imum at a temperature T where £ V {T) — 2£, and be- 
comes nonmetallic, decreases, with further decreasing T. 
The in-plane conductivity a a b remains metallic as long 
as Eq.(17) holds. 

Equations (17) and (18) can be rewritten as the fol- 
lowing relationship between resistivities: 

{PcPabf 12 =P+ (f)pab\ P=4~- (19) 

A fit of the data with Eq. (19) gives p « 0.23 mtt cmpoj. 
Thus, the temperature dependence of p c is determined by 
that of p a b'. 

Note that Eq. (20) predicts a minimum in p c and an up- 
turn at low temperatures. There is a characteristic value 
of p a b = p£o/£> so that when p a b falls below it, p c becomes 
nonmetallic. Since p a b is metallic, and decreases with de- 
creasing temperature, p c is metallic at high temperatures 
and nonmetallic at low temperatures. Specifically, when 
Pab = a T , which is characteristic of the optimally and 
nearly optimally doped cuprates, 

Pc = f3 c + a c T+^. (21) 



The transition from metallic to nonmetallic T depen- 
dence of p c with decreasing temperature is a characteris- 
tic feature of optimally doped and underdoped cuprates. 
In the crystals with maximum T c the onset of supercon- 
ductivity sometimes can mask the upturn in p c so that it 
may look metallic at all temperatures above T c [2l|. This 
is because, according to Eq. (18), p c becomes nonmetal- 
lic when £ v > 2£ and it may not yet reach this value at 
temperatures above T c . 

Using linear extrapolation of p a b — (3+aT with Eq.(20) 
one can predict where the low temperature upturn in p Cl 
normally hidden by the onset of the superconductivity, 
begins. Experimentally, this minimum in p c can be re- 
vealed by the suppression of superconductivity by mag- 
netic field similar to Refs.fll 122. I2I W{ . A possible 
alternative to magnetic field application are the measure- 
ments of the normal state part of tunneling I-V curves 
in mesa junctions, see for example R.efs . [2fiL l27l | . provided 
that that the problem of Joule heating can be adequately 
addressed j2^- More detailed discussion of the applica- 
tion of intrinsic tunneling method to experimental explo- 
ration of the interlayer coherence of single electrons is 
given below in Sec. VII. 

In other types of crystals we may find the situation 
opposite to that in optimally doped Y BaiCu^Oi ~ y . In 
most experimental observations of the normal state resis- 
tivity the temperature range do not extend much above 
the room temperature and, in this range, p c may demon- 
strate the nonmetallic T— dependence. This is typical of 
the Bi 2 Sr 2 C'aCu20g, +x crystals. Optimally doped and 
even slightly overdoped crystals with T c in the range 
80 — 90 K have nonmetallic p c at all temperatures be- 
low room temperature. In contrast, the in-plane resis- 
tivity is metallic and increases with increasing tempera- 
ture. According to Eq.(20) this situation corresponds to 
Pab < P~£o/£,- However, this limit can be reached at suffi- 
ciently high temperature and beyond that p c will demon- 
strate metallic T— dependence. The pronounced mini- 
mum in the T— dependence of p^ in Bi^Sr^CaCu^O^x 
crystals have been reported in|2!| at T ss 750 K. Thus, 
the apparently different temperature dependence of p c in 
optimally doped crystals such as Bi 2 Sr2C aCuiO%+ x and 
Y Ba-iCuzO-j-y may, in fact, be the same, described ap- 
proximately by Eq.(20), except that the minimum of p c 
is either hidden by the onset of superconductivity or lies 
above the temperature range within which the resistivity 
is measured. 

A nontrivial question raised by this analysis is the na- 
ture of the cutoff £ in the ballistic regime, Eq. (17). 
If we assume quasiclassical in-plane transport, namely, 
&ab oc t v , the empirical Eq. (17) gives £ v = q~ 1 o- a b + £ = 
vft v + £. Therefore, the cutoff £ indicates that the in- 
plane phase coherence length does not scale to zero with 
decreasing a a b and t v . I emphasize that this is the high 
temperature cutoff, important only when the decoherence 
time is relatively short. On the other hand, if £ = then 
Cab = q£<p, v c = q£o/£<e and PcPab = const = (g 2 ^j)~\ 
instead of Eq.(19). Without the finite cutoff £ in Eq.(17), 
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the out-of-plane resistivity of incoherent crystals cannot 
have metallic T- dependence at high temperatures. 

The data[2£j that lead to Eq. (17) seem to indicate 
a phenomenon which can be described as a "hard core" 
of the phase-coherent volume. Even when o~ a b —> at 
high temperature, the in-plane phase coherence is re- 
tained within the area of the size ~ £ 2 ~ 20 x 20A 2 . One 
might speculate that such an unusually large value of the 
coherence length at room temperature, accompanied by 
decreasing conductivity, is due to backscattering resulted 
from static disorder. In this sense, the hard core of the 
phase-coherent volume of extended states may have the 
same origin as the phenomenon of Anderson localization. 



IV. 



SCALING APPROACH 



A phcnomenological description of conductivities that 
does not rely on specific microscopic models can be 
developed on the basis of the one-parameter scaling 
hypothesis 0, [ljj. It asserts that the rate of change of 
the conductance of a microscopic block with its size de- 
pends only on the conductance itself and nothing else. By 
extension, the derivative of the conductance of the phase- 
coherent volume dhig v /dln£ v should also be a function 
of g v only, namely: 



din 



9<p 



din L 



(22) 



A nontrivial content of this statement is that the other 
factors affecting conductance, such as the concentration 
of dopants or impurities do not change the functional 
dependence n(g v /g), but may only affect the normaliza- 
tion constant g. It was further suggested inj^,[^ that g 
is universal (~ e 2 /h). This assumption seems to be too 
narrow and cannot hold for all systems. A counter exam- 
ple is a crystal with dopants that change the number of 
carriers, while the coherence length is dominated, let us 
say, by inelastic electron-phonon interaction. Thus, the 
conductance changes with the concentration of dopants, 
while l v remains constant (at a given temperature). For 
Eq. (22) to hold in this case, the normalization constant 
must be concentration dependent. If a crystal is inco- 
herent, the general statement (22) about conductance of 
the phase coherence volume can be translated into ex- 
perimentally verifiable form. 

According to Eq. (9), the anisotropy is a measure of 
the coherence length, £ v = n(o [ry = (<r a b/ 'ct c ) 1//2 ]. On 
the other hand, the in-plane conductivity differs from 
conductance g v by a constant factor, a a b = 9ip/£o (see 
Eq. (5)). Thus, the observable consequence of the one 
parameter scaling hypothesis is that 



d In Ogb 

din T] 



(23) 



In other words, the derivatives d In u a b/d In rj for a family 
of crystals (such as YBa2Cu^Oe+ x , for example) should 



be the same function of o a b and should differ from each 
other only by the a;— dependent normalization constant 
a. For the conductuctivity itself, Eq. (23) translates into 
a two parameter scaling: 



a U 



(24) 



One can speculate that a changes predominantly with 
changing the density of carriers, while fj reflects mostly 
the degree of disorder present in the crystal. The scenario 
presented by Eq. (23) is illustrated in Fig. 1(a). The 
bold segments are schematic representations of n(cr a b) for 
several levels of doping. The range of a a b corresponds to 
"accessible" temperatures T m i n < T < T max . (T max is 
typically 300 — 350 K and T min > T c in superconduct- 
ing crystals). The thin lines indicate the hypothetical 
extensions of the trajectories to "experimentally inacces- 
sible" range of conductivities and temperatures. One 
of the most interesting points here is how the out-of- 
plane conductivity evolves with doping. By definition, 
°~c(t)) = °~ab(v)/ r ] 2 [ see also Eq. (10)], so that 



dln< 



G?ln?7 



= k- 2. 



(25) 



Therefore, both a ab and er c are metallic for k > 2, while 
metallic a a b and nonmetallic a c coexist for < k < 2. 
Both conductivities are nonmetallic for n < 0. Segment 
1 in Fig. 1(a) represents a regime where both conductivi- 
ties a a b and o~ c are metallic at all temperatures T > T m i n , 
because the whole segment is located above the threshold 
k = 2. 

Segment 2 corresponds to a slightly underdoped sys- 
tem. At high temperatures (smaller o~ a b) k > 2 and, 
therefore, a c is metallic and increases with increasing 
t] oc £p\ a c reaches a maximum when K(o~ a b) = 2 and 
decreases with further increasing £ v . Therefore, at this 
level of doping a c has a maximum (p c has a minimum) 
within the accessible T- interval (T m j n < T < T max ). The 
in-plane conductivity remains metallic. 

Segment 3 represents a moderately underdoped sys- 
tem. It is located entirely within the range < k < 2 
and corresponds to metallic o a b and nonmetallic a c for 
all T m i n <i T <C T niax . 

Segment 4 corresponds to strongly underdoped crystals 
with the in-plane conductivity changing from metallic at 
high T (k > 0) to nonmetallic at lower T (k < 0). The 
singularity k = is integrable: 



k(ct) w ±C In 



1/2 



so that o~ab{v) determined by the equation 
r^o din a 



In 



m 



(26) 



(27) 



reaches the maximum value a max at a finite rj = r]i 
and decreases with further increasing r\ (decreasing tem- 
perature). Equations (26) and (27) give cr ab (?]) w 
a max exp{-C 2 ln 2 (ry/?7i)/4}. 
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FIG. 1: (a) Sketch of the trajectories n(a a b) vs. \na a b- Bold 
segments correspond to "experimentally accessible" values of 
conductivities as explained in the text. The thin lines are 
hypothetical extensions of the trajectories drawn under the 
assumption that the segments corresponding to different dop- 
ing levels are parts of a continuous curve k(ct/<t), shifted with 
respect to each other due to different normalization constants 
a. (b) Sketch of the insulating (A) and metallic (C) branches 
of the universal trajectories. Trajectory (B) is the separatrix 
and (D) is the hypothetical "supermetallic" branch. 



Finally, segment 5 corresponds to an insulating crystal 
with nonmctalllic a a b and a c . As shown, all five curves in 
Fig. 1(a) represent the same dependence n{a/a) shifted 
with respect to each other, when plotted against logo - , 
because of the different values of a which is determined 
by the density of carriers n and the level of disorder. The 



reduction of n by doping reduces a, shifting the respec- 
tive segment of n(a /a) to lower values of a and, at the 
same time, ahead along the trajectory in terms of the 
reduced variable a/a. Thus, the experimental data plot- 
ted as a vs anisotropy would present the "snapshots" of 
different stages of evolution of a given system along the 
common trajectory. For example, close to optimal dop- 
ing, p c is metallic in yjSo^Ctt-iO^pH. presumably be- 
cause the minimum in p c is hidden by the onset of super- 
conductivity. The decrease of oxygen content shifts the 
"observable part" of the trajectory so that it crosses the 
threshold n = 2 and the minimum in p c becomes evident 
at T c < T < 300 K. This corresponds to YBa 2 Cu 3 O x 
with x ~ 6.88 and 6.78 which exhibit a minimum in p c 
at T w 150 and 300 K, respectively 30] . At even lower 
oxygen content, the minimum in p c shifts above the room 
temperature and p c appears nonmetallic at all tempera- 
tures below 300 K. 

A single universal trajectory shown in Fig. 1(a) pre- 
dicts that, as temperature decreases, the incoherent crys- 
tals at all doping levels would eventually become insulat- 
ing. First, p c becomes nonmetallic when k < 2 and then, 
at even lower temperature, p a b turns nonmetallic (k < 0). 
Support for this scenario can be found in the experiments 
where superconductivity was suppressed by the magnetic 
field to reveal the underlying normal state. In the sam- 
ples of Bi2Sr2CuO y ^j^ and Bi2Sr2- x La, x CuOe+s^M 
the prominent upturn in p ao has been revealed below T r . 
The same phenomenon takes place in Lai-^Sr^-CuOA |24| 
and Prr>-^Ce^CuOA\2f%. 

The results of the search for the universal trajec- 
tory described by Eq.(24) was reported in RefijU- We 
have found that when the conductivity is plotted vs. 
anisotropy, the appropriately normalized data for several 
crystals of YBozCuzOq+x and Yi^ x Pr x Ba2Cu^OT-s 
with widely different composition fall on the same curve. 

The analysis of the experimental data also suggest that 
one type of trajectory n(a) shown in Fig. 1(a) cannot de- 
scribe all incoherent crystals. Just like the genuinely 2D 
system j32^. the incoherent crystals may undergo metal- 
insulator transition. This scenario is schematically de- 
scribed in Fig. 1(b) by the trajectories A, B and C. The 
trajectory A is the same as in Fig. 1(a) and corresponds 
to the insulating phase. The trajectory C describes the 
metallic branch, and B is the separatrix. The notion of 
insulating and metallic phase refers only to the in-plane 
resistivity. The out-of-plane resistivity is nonmetallic in 
either case as long as k < 2. In the limit n — > the sep- 
aratrix can be approximated as k ~ a\na res /a. Then, 
Eq. (27) gives 

(Jab = Ores exp{-?7~ Q } » a res (l - if a ) (28) 

Taking into account the relationship between the 
anisotropy and coherence length, one can rewrite the Eq. 
(28) as follows: 

Pab = Pres exp j {jf^j j ~ Pres (l + {jf^j ) ( 29 ) 
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FIG. 2: Conductivity of Yi- y PryBa2CuzOj-s single crystals vs. logarithm of anisotropy. The solid lines are guide to the eye 
that correspond to a a b oc lnr; dependence. The inset shows conductivity vs. logarithm of temperature. 



The out-of-plane resistivity diverges in the same temper- 
ature range. From Eq. (28) folows: 



Pab 



, , Pab 

- Pres exp \ [ — 



a/2' 



which gives 



Pab 



1+2/c 
Pres 



[Hpab/presW /a {pab - Pres) 2/a ' 



(30) 



(31) 



Typically, the anisotropy r\ is very large and p a b de- 
scribed by Eq. (28) is practically temperature indepen- 



dent, equal to the residual value p r 



One of 



the samples of Bi 2 Sr 2 CuO y studied in Ref.|23j appears 
to exhibit the properties of the " separatrix crystal" . Its 
in-plane resistivity remains temperature independent be- 
tween 11 K and 0.7 K. The out-of-plane resistivity in 
the same range of temperatures diverges. 

Next, let us consider the metallic branch, the trajec- 
tory C in Fig. 1 (b) . It describes p a b that is metallic at all 
temperatures (k > for all values of a ao ). In incoherent 
crystals the phase-coherent volume is "two-dimensional" 
(contains only one bilayer) and therefore, in the limit 
£ v — ► oo, the conductance should become size indepen- 



dent and, correspondingly, 
k can be approximated as 



. Let us assume that 



At(cr) w a I a. 



(32) 



In a different context a similar asymptotic dependence 
of the logarithmic derivative was discussed in Ref . |33| . 
Integration of Eq. (27) leads to 



ab = a 111(77/170), 



(33) 



where 770 is a constant of integration (77 > 770). Given 
the definition of 77 = {p c / Pab) 1 ^ 2 , this translates into the 
following relationship between resistivities: 



VoPab exp 



P 

Pab 



2a' 



(34) 



Both parameters 770 and p are doping dependent. While 
p a b is metallic at all temperatures, p c reaches minimum 
at pab = P- Depending on where p lies with respect to 
the range of measured p ao , we will see all three types of 
p c (T) dependence discussed above. At sufficiently low 
temperatures the temperature dependence of the coher- 
ence length is a power law, which translates into the log- 
arithmic increase of the conductivity: 



a ab = o-ln.(4,/4)»7o) ~ ln(Xb/T), 



(35) 



while the out-of-plane resistivity p c diverges as the power 
law. 

In Fig. 2 the logarithmic dependence given by Eq. 
(33) is illustrated by the data obtained on several sam- 
ples of Yi-yPr y Ba2Cu^OT-s. This is the same data as 
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p ab ~ ln(l / T) 



FIG. 3: Sketch of three branches of the scaling trajectories: insulating [non-hopping) (A), metallic (C) and separatrix (B). 
Unlike Fig.l, ft(a)is shown vs. a, not In a. 

in Ref. |3ll|. In |3ll | it was presented as lncr vs. In 77 in or- 
der to demonstrate the scaling given by Eq. (24), but 
such presentation obscured the functional dependence 
f(j]) oc In r\. The inset to Fig. 2 shows the conductivity 
of the same crystals vs. logarithm of temperature. It 
is obvious that the conductivity cannot be described by 
the logarithmic temperature dependence. This is because 
the T— dependence of l v in this temperature range has 
not yet settled into power law, but the trajectory k(ct) 
is already well described by the asymptotic dependence 
(32). Apparently, even the n > part of the trajectory 
of the insulating sample YQ ^PrQ ^Ba2Cu^OT^s is well 
described by Eq. (32). 

Let us now turn to the insulating branch, the trajec- 
tory A in Fig. 1(b). Usually it is assumed that in the 
insulating regime the mechanism of conduction is hop- 
ping: a ~ exp{— %£ip}, which leads to k(<t) = ln(cr/<7o). 
There is, however, another possibility. In 2D insulator 
the limiting value of k(ct) may be zero, just as for the 



metallic branch: ft(er) — > as t v — > 00 and a — -> 0. In Fig. 
3 the alternative trajectory of insulating (non-hopping) 
branch (A) that illustrates this scenario is shown. As- 
suming that k{ct) is an analytical function, the simplest 
option is 

K (a) w -— (36) 
Co 

when a — > 0. This leads to 

Pab = po InfaAto )■ (37) 

At low temperatures, when £ v ~ J 1-7 , resistivity will 
acquire the logarithmic T— dependence: 

Pab = po m(V4?7o) - HTo/T). (38) 

This logarithmic T— dependence of the resistivity ob- 
served in Bi2Sr2-xCuOa + s^^ and Y BaiCuiQ7- v \n\ 
has attracted a great deal of attention and has not been 
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explained by any microscopic model. The trajectory A 
in Fig. 3 describes the in-plane resistivity that has metal- 
lic temperature dependence at high temperature, reaches 
minimum when n = and logarithmically diverges at 
low temperature. The out-of-plane resistivity diverges as 
a power of temperature 11] . The trajectories B and C 
in Fig. 3 describe the separatrix and metallic branch, 
respectively, that are the same as in Fig. 1(b). 

The last trajectory (D) shown in Fig. 1(b) describes a 
" supermetallic" branch, such that both p c and p a b are 
metallic at all temperatures (k > 2). If we take the 
asymptotic behavior of K ~ 2 + cro/a, the resulting 



crab = -z- — 
2 \Vo 



2 ' 



(39) 



Here o\ and t\q are constants of integration. Using the 
definition of rj we get 



Pc 



IPab- 



(40) 



While p a b —> 0, p c tends to a finite value. The asymptotic 
value k = 2 corresponds to diffusive in-plane motion of 
carriers. If we take the quasiclassical a a b ~ t v and I 2 ~ 



Dt v , then a a b ~ which corresponds to k = 2. 



V. MAGNETORESISTIVITY 

The effect of the magnetic field on resistivity of the in- 
coherent crystals can be straightforwardly deduced from 
Eqs.(9) and (10) as long as we know what effect the mag- 
netic field has on the coherence length. Hereafter I will 
discuss only the orbital part of the magnetoresistivity. 
Spin-dependent and orbital contributions can be sepa- 
rated due to their different field dependences|34l 133. l3q . 
Normally, the application of the magnetic field perpen- 
dicular to the conducting layers decreases 4, 131 . The 
first obvious observation that follows from Eq.(9) is that 
magnetoanisotropy is always negative: 



Pc/ Pab 



L 



< 0. 



(41) 



Magnetoanisotropy can be measured directly and, in 
fact, more accurately than the separate magnetoresis- 
tivities Ap a b and Ap c by the six-point methodj^ E3- 
Even when both coherence lengths (in-plane and out- 
of-plane) change with temperature, the Eq.(41) holds 
as long as the the coherence length in the direction 
of the applied field does not change by the field, see 
Eq.(6). Indeed, in Refs.jU 13 the negative magne- 
toanisotropy was observed in all samples of Bi 2 Sr2CuO y 
and Y Ba2Cu^Oj-s- 

Since both conductivities, a a b and a c are determined 
by the in-plane coherence length, Eqs.(23) and (25), there 
is a direct relationship between the type of conduction 
(metallic or nonmetallic) and the sign and magnitude of 



the magnetoresistivities (MR): 

Apab_ M v Ap c _ A£ v 

— — ~ K ~p — ' ~ — — \ K Z /p — i 

Pab "ip Pc *-tp 



(42) 



For < k < 2, Apab is positive and Ap c negative. 
In other parts of the trajectories both magnetoresistivi- 
ties are either positive or negative. Since the tempera- 
ture dependence of the resistivities is determined by the 
T— dependence of £ v , there is a relationship between the 
temperature coefficient of the resistivity and the sign of 
magnetoresistivity: 



dpab,c 

dH 



) £>Pab 2 c 

dT 



(43) 



where Q = (de ifi /dH)/(dt lp /dT) > 0. As long as the 
magnetic field reduces the coherence length, its effect 
is equivalent to the increase in temperature. Thus, the 
coexistence of metallic p a b and nonmetallic p c in many 
cuprates translates into opposite signs of magnetore- 
sistivities. The correlation between MR and the sign 
and magnitude of dp/dT has been well documented in 
literaturejH |H H IS H| . 

There is an elegant and conceptually straightforward 
way to verify experimentally the relationship between the 
anisotropy and coherence length, Eq. (9), using the well 
established theory of magnetoresistance due to quantum 
interference. Since the coherent trajectories are confined 
to a single bilayer (and therefore two-dimensional) the 
field dependence of magnetoresistance is given by 



A Pab ,c^^^HB/B v ) 



(44) 



This follows from Eqs.(42) and (43) and At v ~ 
-ln(B/B v )\j3, E|. Here B > B v = hc/AeDr^ = 
(j) /4:nDT v and 4> = nch/e ~ 2 x 10~ 7 G cm 2 is the quan- 
tum of flux. The coherence length in the diffusive regime 



is defined as 



(x 2 ) = (l/2)(r 2 ) = 2Dt v . This follows 



from the probability to find a particle at a distance r from 
the starting point: dW(r,t) cx exp{— r 2 / ADt]dr. Thus, 
Bp = 4 i o/2'k£ 2 . According to Eq. (9), in incoherent 
crystals a product of seemingly totally unrelated quanti- 
ties B v p c / pab should remain constant, even though both 
the anisotropy of resistivity and B v strongly change with 
temperature. Moreover, a combination such as 



£ - ( ^° ° c 
V,C \2ttB v Oal 



1/2 



(45) 



determines the out-of-plane phase coherence length and 
should yield a value of the order of interatomic (inter- 
layer) distance. The logarithmic field dependence of 
MR has been observed in the normal state of layered 
cuprates[!3liUE3- At sufficiently low temperatures the 
value of B v can be reliably established by fitting MR 
data to the theory of weak localization. This type of 
experiment requires both the anisotropy and MR to be 
measured on the same crystal. 
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FIG. 4: Temperature dependence of the anisotropy of two underdoped Y Ba-zCu^Ox single crystals. The anisotropy of 
Y Ba2CuiO6.se crystal increases monotonically in the temperature range shown. The more underdoped Y Ba^CuzOs.^ sample 
has non-monotonic temperature dependence of the anisotropy. The value of p c /p a b ~ 1.9 x 10 3 at T = 75 K is indicated. The 
dashed curve is a hypothetical extrapolation of the high temperature behavior of the anisotropy. It shows what the values of 
the anisotropy would have been, had the out-of-plane coherence length remained constant, equal to 6A, see discussion in Sec. 
VII B. 



In Ref . |34j the magnetoresistivity data obtained on 
YBa 2 Cu30e.25 crystal at T = 75 K indicate that B v w 
0.3 T. The temperature dependence of the anisotropy of 
this sample is shown in Fig. 4 and at T = 75 K it is 
about 1900. Therefore, Eq.(45) yields 



l v , c = 7AA. 



(46) 



Note that in Y Ba^CuaOe^s crystal the anisotropy de- 
crests below w 130 K. This underdoped crystal exhibits 
the properties of the semicoherent system, discussed be- 
low in Sec. VII, and at 75 K £<* c has already increased 
beyond the minimum value £q = 6 A. 



A. Violation of Kohler's rule 

The origin of the Kohler's rule can be illustrated on 
Eq.(42). In the weak field regime 



A4 



(47) 



where oj c is the cyclotron frequency. Thus, according to 
Eq.(42): 



Ap ab 

Pab 



k(uj c t v ) 2 



In Fermi liquids k = 1 and a plot of MR vs. B 2 t 2 pro- 
duces a straight line with T— independent slope pi). In 
incoherent crystals even this modified Kohler's relation- 
ship does not hold because the trajectory k is not a con- 
stant, Fig. 1. Even if we would somehow manage to 
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reliably determine t v , the slope of the plot Ap a b/p a b vs. 
{uj c t v ) 2 will be proportional to k and change with tem- 
perature. 

An important consequence of Eqs. (41) and (42) is that 
the trajectory k{T) or k{<t) can be extracted from MR 
data if the magnetoanisotropy is measured along with 
the components of MR. Namely from Eqs. (41) and (42) 
follows: 



2 A In p a b 
Ala(p c /p ab ) 



(48) 



Thus, the trajectories of the type shown in Fig. 1 or Fig. 
3 can be directly obtained by measuring both compo- 
nents of MR similar to how it was done in Ref.JsIj]. The 
applicability of Eq. (48) does not depend on whether the 
field dependence of MR is quadratic or not. 



VI. EFFECT OF ELEMENTAL 
SUBSTITUTIONS AND RADIATION-INDUCED 
DISORDER 

The analysis given in the previous section can be 
readily extended to another property of layered crys- 
tals: the response of the resistivities to elemental sub- 
stitutions, especially those that replace Cu in the CuOi 
planes. The impurities reduce the in-plane coherence 
length {dt v /dx < 0) by decreasing the elastic mean free 
path and the diffusion coefficient (x- is the concentra- 
tion of substitutions). Therefore, the anisotropy of the 
incoherent crystals always decreases with concentration 
of such impurities, while the resistivities change accord- 
ing to the position of a given crystal on the trajectory, 
Fig.l(a). 



d In p ab 
dx 



d In l v d In p c 



():> 



dx 



-(«-2)^. 



(49) 



This conclusion can be illustrated on the examples of 
Bi 2 Sr ? Ca(Cui X M X ) ? Q H+V (M = Zn, Mn, Fe, Co, 
and iVa)[42l l43j. Undoped BiiSriCaCuiO^+y crystals 
exhibit metallic p a b{T) and nonmctallic p c {T), so that 
below room temperature < k < 2 and, as the result, 
dpab/dx > while dp r /dx < 0. All elemental substi- 
tutions examined in [42L |43j lead to decreasing of p c and 
increasing of p ab - 

Sometimes, the increase of the c-axis conductivity in 
response to a perturbation has been interpreted in lit- 
erature as a crossover to coherent transport in the c- 
direction. We see that this is not necessarily the case. 
The anisotropy and p c may decrease due to the reduc- 
tion of the in-plane phase coherence length, even when 
the c-axis coherence length remains fixed. 

Radiation-induced disorder in the CuOi planes may 
have the same effect on the resistivity as the elemental 
substitutions. The metallic in-plane resistivity of non- 
irradiated crystals and films increases with the amount 
of absorbed radiation and even changes its temperature 



dependence 44] revealing a minimum in p a b at low tem- 
perature. Unfortunately, I am not aware of any report 
on the effect of irradiation on the out-of-plane resistivity. 

Unlike magnetoresistivity, the response to substitu- 
tions for Cu on the planes is very dramatic. Small, 
of the order of one percent, doping drastically affects 
the magnitude and temperature dependence of the re- 
sistivities. Therefore, the differential relation given by 
Eq.(49) can only give a qualitative idea of the changes 
introduced by this type of doping. For example, in 
Bi2Sr2Ca(Cui- x Zn x )20s+ v crystals just 1% of Zn dop- 
ing reduces the anisotropy by an order of magnitude. The 
change in concentration of such dopants may not only 
drive the crystal along the given scaling trajectory, but 
can cause the metal-insulator transition. 

One can ask why a very small concentration of defects 
in the CuOi planes so drastically alters the resistivity 
and the in-plane coherence length, while the response 
to other substitutions, for example to Pr substitution 
for Y, is much more gradual. One scenario is that the 
defects on the planes violate the reflection symmetry of 
the CuOi bilayers and cause hybridization of the even 
and odd subbands pM UK I47 L l48j . If the relaxation rates 
of these subbands are very different, this hybridization 
may lead to drastic reduction of the decohcrence time. 



VII. SEMICOHERENT CRYSTALS 

Between two extremes, conventional Fermi liquids with 
temperature independent anisotropy on one hand and 
incoherent layered conductors, exhibiting the strongest 
possible rate of change of the anisotropy on the other, 
there is an intermediate class of layered crystals. In these 
systems the out-of-plane coherence length changes with 
temperature, but not at the same rate as the in-plane 
coherence length. As the result, in such "semicoherent 
crystals", the anisotropy changes with temperature, but 
not as strongly, and not necessarily monotonically as in 
incoherent layered systems. 

The first simplest scenario of semicoherence is that 
£ip,ab ^> £<p,c and t v , c is so short that the Ioffe-Regel 
limit still affects its temperature dependence: l VtC — 
^o + v f,cT v , while £ Vi ab = VF,abi~ v . This empirical interpo- 
lation form of £ VtC takes into account that the coherence 
length cannot be arbitrarily small and saturates at high 
temperature at a finite limit. The anisotropy then is 
given by 



0_ab_ 

Or. 



V 2 T 2 
v F,ab T v 



{t + VF^Ttp) 2 



(50) 



The anisotropy monotonically increases with decreasing 
temperature and saturates at the level determined by the 
anisotropy of the three-dimensional Fermi surface. If we 



take the conventional form of metallic p a b 
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FIG. 5: Anisotropy vs. in-plane resistivity for several Nd,2- 
straight lines are power law fits to the data. 



c Ce x CuOA- v crystals. The exponent /3 is defined by Eq.(52). The 



out-of-plane resistivity is given by: 



A. Nd,2-xCe x Cu04- y crystals 



Vmax 



Pab 



{Pab/P+If 



(51) 



It is nonmetallic at high temperatures, reaches maxi- 
mum when p a b = p and decreases at lower temperatures. 
There are several types of layered crystals that exhibit 
the temperature dependence of the resistivities qualita- 
tively similar to that described by Eq. (51). In Ref.|49| 
the coherence-incoherence transition was examined with 
angle-resolved photoemission spectroscopy and electronic 
transport measurements. In (Bi 5 Pbo^)2Ba 3 Co20 y 
and NaCo20i crystals the crossover to 3D coherence 
is accompanied my a maximum in p c (T), while p a b(T) 
remains metallic and anisotropy increases with decreas- 
ing temperature and saturates at a large value. A simi- 
lar behavior of p c {T) and p a b(T) characterizes S^RuO^ 
crystals 50] . Note that the maximum in p c in semico- 
herent crystals is in stark contrast to the minimum in 
p c in incoherent crystals, Eqs.(20),(21) or (34). One can 
suggest that in some crystals it might be possible to ob- 
serve the out-of-plane resistivity with two extrema: the 
minimum at higher temperatures where the crystal is in- 
coherent and the maximum at lower temperature that 
indicates a crossover to 3D coherence. 



Another type of scmicohercnt systems that presents 
a greater challenge to interpretation is electron-doped 
Nd 2 - x Ce x Cu04- y single crystals. These layered crys- 
tals have large anisotropy ~ 10 4 — 10 5 even though 
both resistivities p a b and p c have metallic T— dependence 
above T c . The temperature dependence of the in-plane 
resistivity is quadratic as in conventional Fermi liquid 
type metals. In Ref.jHlJ the resitivity of several crystals 
with nominal Ce concentration of x = 0.08 and 0.29 were 
presented. The lower and higher Ce concentration cor- 
responds to underdoped and overdoped regimes respec- 
tively. However, because of uncontrolled oxygen content 
the resistivities vary even within each group with the 
same Ce content. In Fig. 5 the data from Ref.[5l| are 
presented in order to identify the relationship between 
anisotropy and in-plane resistivity. For all samples it can 
be well approximated by: 



Pab 



^Pab 



APab 



(52) 



The 3D coherent transport corresponds to the exponent 
[3 = 0, while the incoherent transport is characterized 
by (3 = 2 (Graf, Rainer, Sauls relationship 8] , Eq. (11)). 
If we take p ab ~ r" 1 - £~] Eqs. (6) and (52) are 
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equivalent to the following relationship between the in- 
plane and out-of-planc coherence lengths: 

VcOc££ /a . (53) 

Note that in Fig. 5 the exponent j3 systematically de- 
creases with overdoping (higher Ce concentration). 



B. Nonmonotonic temperature dependence of 
anisotropy 

The third scenario of gradually developing coherence 
in the c— direction differs from the previous two in the 
way the anisotropy changes with temperature. Both Eqs. 
(50) and (52) describe monotonically increasing with de- 
creasing temperature anisotropy and metallic in-plane re- 
sistivity. In literature one can find several examples of 
underdoped crystals in which p a t is non-metallic at low 
temperatures and anisotropy changes nonmonotonically, 
namely, exhibits a maximum at finite temperature. 

One example of such behavior is shown in Fig. 4. The 
less underdoped Y Ba2Cu3OQ.se crystal has monotoni- 
cally increasing anisotropy. The anisotropy of stronger 
underdoped Y BaiCu?,0§.i : ~ s crystal reaches maximum at 
T m w 130 K and decreases rapidly at lower tempera- 
tures. The temperature dependence of the resistivities of 
these samples is shown in Ref.j23]- The Y Ba2Cu?,0§.z§ 
crystal has metallic p a b, while Y BaiCu^Oe.ih has non- 
metallic p a b at T < T m . 

According to Eq.(6), the decrease of the anisotropy at 
lower temperatures indicates that below T m the out-of- 
plane coherence length increases more rapidly than the 
in-plane coherence length. Presumably, if the crystal re- 
mains in the normal state at sufficiently low temperature, 
or superconductivity is suppressed by the magnetic field, 
the anisotropy of crystals like Y Ba^Cu-iO^.i^ will even- 
tually stabilize at the level determined by the anisotropy 
of the 3D wave functions of the localized electrons, which 
is smaller than the maximum value of the anisotropy at 
intermediate temperatures. 

There are several points worth noting. There is a 
widely held opinion, expressed in literature, that the ten- 
dency to establish the out-of-plane coherence increases 
only with increasing number density of charge carriers, 
while the underdoped crystals are incoherent and be- 
come more so with decreasing number of charge carri- 
ers. Especially in Y Ba2Cu3,O y , the disruption of the 
CuO chains is thought to be the main reason for increas- 
ing anisotropy. This is obviously not the case, and the 
systematic stu dy o f the evolution of the anisotropy with 
oxygen content [3 Ij shows that the anisotropy at a given 
temperature reaches maximum at an intermediate oxy- 
gen content. 

Apparently, in cuprates there is a mechanism of in- 
terlayer decoherence, of still unknown nature, which can 
be overcome by the phase transition into the supercon- 
ducting state. This takes place in optimally and slightly 



underdoped cuprates. If, however, the onset of supercon- 
ductivity is somehow prevented, one can expect that at 
sufficiently low temperatures 3D coherence will still be 
established. As in the case of YBa2CusOe.25, Fig. 4, 
this will manifest itself as the decrease in anisotropy be- 
low a certain temperature. How low is this temperature 
is a measure of strength of the decoherence mechanism. 
Then, it is obvious that the mechanism of decoherence 
is stronger in YBa2Cu3.Os.3Q and other YBCO crystals 
with higher oxygen content than in Y Ba2CusOe.25- This 
observation has to be contrasted with the behavior of 
Pr doped, full y o xygenated, Yi^ x Pr x Ba2Cu30rs sin- 
gle crystals |l8l l3l| . Even the very strongly underdoped 
YoA7Pro.53Ba2Cu 3 C>7^s and insulating PrBa2CusOrs 
remain incoherent, with increasing anisotropy, at tem- 
perature as low as 1.9 K. 

A conclusion that one can make from this is that 
in YBCO and its derivatives the CuO chain layer 
is responsible for the interlayer decoherence. In 
Yi^ x Pr x Ba2Cu'iO'j^s crystals the chain layer is fully 
formed and the mechanism of decoherence is strongest. 
In Y Ba2Cu^Oe.25 this mechanism is weakened to such 
extend that it allows the process of establishing 3D co- 
herence to begin at temperature as high as 130 K . 

The central role of the chain layer in YBCO is also 
corroborated by the value of £ VtC w 7.4 A, Eq.(46). At 
this temperature (75 K) the anisotropy p c /p a b ~ 1900 
and l VtC is already increasing beyond its minimum value 
£q « 6A at higher temperatures. If ^ VjC were still 
equal 6A at T — 75 K the anisotropy at this temper- 
ature would have been 1900 x (7.4/6) 2 = 2900. This is 
close to the value of the anisotropy one gets extrapolat- 
ing the high temperature behavior of the anisotropy of 
YBa2Cu3.OQ.2b as shown by the dashed line in Fig. 4. 

Thus, the loss of phase coherence of the charge carriers 
takes place during their transfer from the CUO2 planes 
to the charge reservoir. Otherwise, if the loss of coher- 
ence was taking place during the charge carriers trans- 
fer between the neighboring CUO2 bilayers, the value of 
£0 would be equal to the distance between the bilayers 
(12A) (see Section III and comparison with the solvable 
modelQ). 



C. Intrinsic tunneling spectroscopy in 
semicoherent crystals 

The intrinsic tunneling spectroscopy in small 
crystals!!! and mesa-structures^ Km Ei.EEiE^ 
may provide a powerful tool for investigation of the 
process of establishing the interlayer coherence of single 
charge carriers. The most studied crystals such as 
Bi 2 Sr2CaCu20 8 +x (Bi - 2212) and BiiS^CmOs+s 
(Bi — 2201) are found to have current-voltage (TV) 
characteristics similar to that of a stack of Josephson 
junctions. These intrinsic junctions are formed by the 
charge reservoir layers sandwiched between supercon- 
ducting CUO2 units (a superconducting unit is either a 
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single CuOi layer like in Bi — 2201 or a CuOi bi-layer 
in Bi — 2212 and YBCO). Below critical temperature 
the I-V characteristics consist of multiple hysteretic 
branches. Each branch corresponds to switching of 
one intrinsic junction from the superconducting to the 
normal state. 

In mesa structures the number of charge reservoir lay- 
ers can be made small (10 — 20) and, correspondingly, 
the number of branches is easily countable. For example, 
in Ref.|26j| the number of branches remains the same in 
the temperature interval from 70 to 4.2 K (Bi — 2212). 
This is understandable if the single quasiparticles remain 
incoherent, so that the Cooper pairs from each individ- 
ual CuOi unit form a Josephson junction across every 
charge reservoir in the mesa. The normal state c— axis 
resistance of that sample increases with decreasing tem- 
perature, corroborating the assertion that this crystal is 
incoherent. 

What would change in this picture if the crystal was 
semicoherent? If the out-of-plane phase coherence length 
of single quasiparticles £ VjC is already greater than its 
minimum value £q at temperatures above T c , but at T c 
is still shorter than the height of the mesa, the num- 
ber and the properties of the intrinsic junctions will be 
different from those formed in incoherent crystal. To il- 
lustrate that, let us say that just above T c , in the normal 
state, a number of neighboring CuOi units have formed 
"coherent clusters". This means that the single charge 
quasiparticles in these clusters consist of the orbitals of 
several (two or more) neighboring CuOi units, as op- 
posed to incoherent crystals where the wavefunctions of 
the quasiparticles consist of the orbitals of one CuOi unit 
only. The rest of the CuOi units in the mesa remain in- 
coherent. 

Let N be the number of CuOi units in the mesa. And 
let rii be the temperature dependent number of double 
coherent clusters (712 < N/2). Then N — 2ni is the num- 
ber of single incoherent units. For simplicity we have 
neglected the probability of formation of coherent clus- 
ters consisting of three and more CuOi units. 

Below T c the quasiparticles form Cooper pairs in both 
coherent clusters and single units. When the supercon- 
ducting order parameter is extended throughout the crys- 
tal, the intrinsic junctions with the lowest critical cur- 
rent J c are the ones between two neighboring coherent 
clusters, or between a single unit and coherent cluster, 
or between two single units. The total number of these 
junctions is N — n%(T). In other words, a charge reser- 
voir "hidden" inside every double cluster does not form a 
weak junction. Thus, the number of hysteretic branches 
in I-V characteristic N\, = N — niiT) and, therefore, can 
be noticeably smaller than the number of unit cells along 
the c-axis. Even more telling, the number of branches 
will be temperature dependent, decreasing with decreas- 
ing temperature. 

Generally, we can introduce nk as the average number 
of coherent clusters consisting of k CuOi units so that 
knk = N. The number of charge reservoir lay- 



ers hidden inside these clusters is 53fc=i(^ — l) n fe- Then 
the number of low— J c hysteretic junctions and, corre- 
spondingly, the number of low— J c hysteretic branches 
Nb = rifc - the total number of coherent clusters. 

In incoherent crystals Nb = N, since the only clusters are 
single units. As the coherence of single charge quasipar- 
ticles sets in, the number of coherent clusters along the 
height of the mesa decreases and in the limit when (.^^ 
exceeds the height of the mesa, the mesa itself becomes 
one coherent cluster so that only one branch of I — V 
characteristic, V = 0, remains as long as I < J c . 

Obviously, this a very simplistic description of a com- 
plex process of breaking up of coherence along the long 
chain of CuOi units. For example, one can expect that 
once the lowest— J c junctions are broken, the increasing 
current will start breaking coherence of the coherent clus- 
ters. This may result in another sequence of hysteretic 
branches. Joule heating will complicate further this pic- 
ture. My point, however, is that in the mesas made of 
semicoherent crystals the pattern of hysteretic branches 
and its evolution with temperature might be distinctly 
different from that found in Bi - 2201 or Bi - 2212. 

Semicoherent crystals for this type of experi- 
ment, besides those already mentioned above, may 
include RBaiCu^ Ofu - x (R = Tm. Lu) |57| and 
Ndi- x Ce x Cu04- y 58]. A characteristic feature that in- 
dicates semicoherence is the maximum in p c or anisotropy 
at T = 100 - 150 K. 



VIII. CONCLUSIONS AND SUGGESTIONS 

This paper presents a phenomenological approach to 
understanding the normal state transport properties of 
incoherent and semicoherent layered crystals. It is based 
on the fundamental relationship between the resistive 
anisotropy and the phase coherence lengths derived in 
Sec. II. This relationship is especially useful when one of 
the coherence lengths is fixed, temperature independent 
distance as is the case in some layered crystals such as su- 
perconducting cuprates and many others. In Sec. Ill we 
have shown that the results obtained by our approach are 
equivalent to those obtained for a solvable microscopic 
model Sj. Application of this approach to hopping con- 
duction in layered crystals and to crystals like optimally 
doped cuprates show a good agreement with experimen- 
tal data. 

In incoherent crystals the resistive anisotropy is a mea- 
sure of the in-plane coherence length and, therefore, such 
systems allow an effective application of the scaling the- 
ory as demonstrated in Sec. IV. One of the important 
consequences of the scaling approach is the idea of the 
universal trajectories. We suggest that the families of 
cuprates such as Y BaiCu^O^+ y with different oxygen 
content or Y\- x Pr x BaiCu-z07-s with different Pr con- 
centration are described by the same dependence of con- 
ductivity vs. anisotropy, Eqs. (23), (24). This type of 
scaling has been found in Ref|3l|. 
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Moreover, the scaling approach allows to predict two 
important asymptotic branches of the in-plane conduc- 
tivity. Using simple analytical form of the scaling trajec- 
tory we have shown that the metallic branch of incoher- 
ent crystal can be described by logarithmically increasing 
conductivity: o~ a b oc ln(^ y ). This asymptotic behavior is 
illustrated in Fig. 2. For the insulating branch, using a 
simple analytical expansion of the scaling trajectory, we 
obtain the resistivity that is increasing logarithmically: 
p ao oc hx(£u}). At sufficiently low temperatures, when 
£ v ~ T~ 7 , this crosses over into the famous logarithmic 
temperature dependence: p a b oc ln(l/T). 

In Sec. V the magnetoresistivity of incoherent crys- 
tals is discussed. One of the most interesting conclusions 
is that at temperatures where MR due to quantum in- 
terference can be observed one can also determine the 
out-of-plane phase coherence length, Eq.(45). In Sec. 
VI we consider the effect of elemental substitutions and 
radiation induced disorder. The treatment used in the 
previous sections also allows to explain sometimes puz- 
zling conflicting responses of in-plane and out-of-plane 
resistivities to the perturbations such as impurities and 
imperfections. 

In Sec. VII we discussed the semicoherent crystals in 
which both in- and out-of-plane coherence lengths change 
with temperature, but at different rates. Three scenarios 
are discussed; all of them, it seems, can be found in either 
overdoped or underdoped cuprates. 

A. Possible future experiments 

Typically in literature the problem of the normal state 
of the cuprates is stated in terms of the anomalies of 
the out-of-plane conductivity. We see that the proper- 
ties of o~ c are entirely determined by the fact that the 
out-of-plane coherence length has a fixed value and does 
not change with temperature, Eq. (10). Then, the real 
unresolved question about the anomalous properties of 
the normal state in cuprates is the nature of the strong 
mechanism of decoherence of the charge carriers over very 
short distance (confinement). While this question is im- 
portant in and of itself, it is also a key question which 
is necessary to address in order to understand the super- 
conductivity in cuprates, which is likely to be the phase 
transition that allows to overcome the confinement. 

One possible scenario of decoherence was outlined in 
Ref . [5^| . where the confinement was attributed to the 
dephasing of the tunneling charge carriers resulted from 
the interaction with charge fluctuations. While this effect 
appears to be too weak to account for the phenomenon 
of confinement, it raises a question: how strong is the 
dephasing mechanism? For example, as was discussed 



above, Pr doped YBCO crystals remain incoherent even 
at T k 2 K . However, if the confinement is the result of 
a conventional process of finite strength, one can reason- 
ably assume that eventually, at sufficiently low tempera- 
ture, even these systems will begin to form 3D coherence 
in the normal state. 

As was discussed in Sec. VII B the decoherence mech- 
anism appears to be determined by the properties of the 
charge reservoir and is weaker in both overdoped and 
underdoped cuprates. In underdoped cuprates the tran- 
sition to 3D coherence manifests itself as the maximum of 
anisotropy at finite temperature which has been observed 
in YBCO, Fig. 4, and in Nd 2 - x Ce x Cu0 4 - y crystals Hj. 
However, the strongest confinement in YBCO is asso- 
ciated with fully oxygenated CuO chain layer. In my 
view, it would be very important to find out if, indeed, 
the fully oxygenated derivatives of YBCO, in which su- 
perconductivity is suppressed, will exhibit the maximum 
in anisotropy in the normal state. This can be done with 
Pr doped YBCO or Y Ba2{Cui- x M x )^OT-s, or a com- 
bination of both types of doping optimized with the goal 
to minimize T c , but retain reasonably high conductivity. 

Another important experiment would be to verify the 
Eq.(45). Since the theory of MR due to quantum inter- 
ference is well established, this relationship is a direct 
consequence of our main result, Eqs. (6) and (9). At this 
point it has been applied at only one temperature and the 
result looks reasonable, Eq. (46). If it can be substan- 
tiated for different types of incoherent crystals that the 
length determined by Eq. (45) is indeed of the order of in- 
terlayer distance and remains relatively T— independent 
over a wide range of temperatures, it will be an exper- 
imental proof of the fundamental relationship between 
the resistive anisotropy and coherence lengths. 

There remains also the question of the universal tra- 
jectories, Eq. (23). In Ref.|3l| the two parameter scaling 
dependences 17(77), Eq. (24), were obtained by measuring 
the resistivities p c and p a b- However, the magnetoresis- 
tivity measurements allow to go further and obtain the 
one parameter differential trajectories k(o~) directly ac- 
cording to Eq. (48). 

Last but not least is the possibility to observe the 
establishing of interlayer coherence with the help of 
intrinsic tunneling spectroscopy. To my knowledge, no 
experiments on the crystals in which p c (T) exhibits a 
crossover from incoherent to coherent transport have 
been conducted so far. 
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